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j/3 ■ In this paper, we reconstruct cosmological models in the framework of f{R,T) gravity, 

o: 

^~>, the dust fluid reproduces ACDM, phantom-non-phantom era and the phantom cosmology. 

^. 
P> f Further, we reconstruct different cosmological models including, Chaplygin gas, scalar field 

with some specific forms of f{R,T). Our numerical simulation for Hubble parameter shows 

^ ■ good agreement with the BAO observational data for low redshifts z < 2. 

^ ; PACS numbers: 04.20.Cv, 04.50.Kd, 98.80.Jk, 98.80.Bp 

in 

o 



where R is the Ricci scalar and T is the trace of the stress-energy tensor. We show that 
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INTRODUCTION 



From the cosmological observational data 



J], it is now well established that the present ob- 



servable Universe is undergoing an accelerating expansion. While the source driving this cosmic 
acceleration is known as 'dark energy' its origin has not been well understood yet due to absence 
of a consistent theory of quantum gravity. This acceleration is driven by the negative pressure of 
the dark energy. The 'cosmological constant' is the most simple and natural candidate for explain- 
ing cosmic acceleration but it faces serious problems of fine-tuning and large mismatch between 

MM 

theory and observations [5-7]. Hence there has been significant development in the construction 
of dark energy models by modifying the geometrical part of the Einstein-Hilbert action. This 
phenomenological approach is called as the Modified Gravity which can successfully explain the 
rotation curves of galaxies, the motion of galaxy clusters, the Bullet Cluster, and cosmological 
observations without the use of dark matter or Einstein's cosmological constant |8l4l9l|. The f{R) 



theories can produce cosmic infiation, mimic behavior of dark matter and current cosmic acceler- 
ation, being compatible with the observational data 
gravity and its cosmological implications ). 



2Ch23|| (also see a recent review |2j] on f{R) 



In a recent paper 



25|, the authors considered a generalized gravity model f{R,T), with T 



being the trace of stress-energy tensor, manifesting a coupling between matter and geometry. By 
choosing different functional forms of /, they solved the dynamical equations of astrophysical and 
cosmological interest. In the present work, we study the same model by taking different kinds of 
the energy sources. Reconstruction of the cosmological models for this theory have been ivestigated 
by several authors [3Ch34 | . 

The paper is organized as follows. In section II we present the general action and the equation 
of the motion for f(R,T). In section III we reconstruct the cosmological models from the dust 
fiuid. In section IV we generalize the dust fluid models to the general fluid with EoS p = oop. 
Section V is devoted to discussion and conclusions. 



II. FIELD EQUATIONS IN f{R,T) GRAVITY 



The action of /(i?, T) gravity is of the form 25 1 



S=r^ I d^Xy/^f{R,T)d^x+ f L^./^d*x, 



(1) 



where f{R, T) is an arbitrary function of the scalar curvature R = i?^ and the trace T = T^ of the 
energy-momentum tensor T^y. We define the Lagrangian density for matter field Lm by 

The equation of motion (EOM) is obtained by varying the action ([T|) with respect to g^'^ 25|] 



fuRfM,^ - ^RgfMu + ig^,uO - V^Vy)fR = S-kT^u - frT^u - h^fiv, (3) 

where fR{R,T) = ^^ and fT{R,T) = -^^^ ' . V^ is the operator for covariant derivative and 
box operator (or d' Alembert operator) D is defined via 

1 „ . , — ,„,„ . „ „,fl5r„ 



Performing a contraction of indices in ([3]) , we obtain 

R/r + SafR - 2/ = SvrT - TfT - Qfr- (4) 

Here = g^'^Q^u- We will use ([4]) to reconstruct f{R,T) for different kinds of matter. 

III. RECONSTRUCTION OF fiR,T) USING DUST 

We consider the metric of a flat Friedmann-Robertson- Walker (FRW) spacetime 

g^y = diag{l, -a\t), -a\t), -a\t)). (5) 

Following the conservation of the energy momentum tensor T-^ = for metric ([5]) we obtain 

t = -3HT. (6) 

The 00— component of ([3]) reads 

/rQr - ij) + 1/ = ^T{8Tr + /t) + R^fRR + 9H^T^fTT - QHTR/rt - 3r(| - 5H^)fT. (7) 

In ([7|), if we know the functional form of the Hubble parameter H as a function of the two in- 
dependent variables R and T, then reconstruction of the model f{R,T) is straightforward. But 
always, we have a physical intuition about the form of H{R,T). We clarify it here more. In a 
FRW model, the form of the Ricci scalar is i? = 6{H -\- 2H^). Any cosmological era has a specific 
matter fields, radiation, phantom field or a mixed of them. For example in ACDM model, we have 
H"^ = Hl + ^^-2|2 — . It is easy to see that we can write R in terms oi H as H"^ = ^ - 3H^. Thus 
usually, we can write H^ = H^{R). Thus (11) can be integrated easily to give the functional form 
f{R,T). Also, since f = -3H{R)T, f = -3r(| - 5iJ2) thus all terms of the © are functions 

of R, T. For term as R we can write it as i? = -^ = — — jjj . Mathematically, Q is a second 

order partial differential equation for f{R,T). 



A. Cosmological Implications 

We now discuss the solutions of ^ relevant in cosmological context. In the coming sections 
we will show that, any cosmological epoch (radiation, matter, dark energy dominated eras) can be 
constructed in a model of f{R,T) only with dust fluid as the source. We mention here that the 
form of the action f(R,T) is not unique. For this reason we can continue our investigation of the 



cosmological reconstruction for another model with f{R,T) = fi{R) + f2{T) 25|. When we set 
/2(T) = 0, we recover the f{R) theory which has been discussed in ]2T\. Adopting the technique 
of the reconstruction of f{R, T) models, we have the following different models reconstructed from 
the dust fluid. 

B. ACDM model 

In Einstein gravity, the Hubble parameter for a flat FRW model with real matter field describes 
the ACDM model by H^ = H^ + ^poa'^, in units k^ = SvrG c = 1. Since 



H' = §-3Hl (8) 

R = -9[^-3Hif\ (9) 

R = ^iR-12Hi){R-9Hi), (10) 



then the field equation ([7]) converts to the following form 



^lR-l(R-12Hi)(^R-9Hi 



fR + lf = InSTT + fr) + 8l(f - SHlYfRji (11) 



+9(1 - miy^JTT + 54(1 - mlfrfRT - 3r(| - bH'')fT. 

The solutions of ()lll) can be obtained under two special cases as: 

• f{R,T) = F{R) and F{R) = AR + B where A, B and T are constants. Equation (15) gives 
B = 327rr and A = 0. 



f{R,T) = G{T) and G{T) = CT + D where C, D and R are constants. Equation (15) gives 

L6/37r 
R+9H^ 



C = iJt/lL^ and D = 0. 



C. f{R,T) reproducing the system with phantom and non-phantom matter 

In Einstein gravity, a Universe filled with a mixture of the phantom and non-phantom com- 
ponents obey from the following expression for H, H^ = ^{pqa~^^ + Ppof^^)- Here the set of the 



variables pq (energy density of non-phantom matter) pp, (energy density of phantom matter) ci are 
the parameters of the model. In the early Universe when the scale factor was small, the first term 
in Pqtt"'^^ dominates and it behaves as the Universe described by the Einstein gravity with a matter 
whose EoS parameter is w = —1 + ci/3 > —1 which means that the matter field is non-phantom 
like. But for present era we have w = —1 — ci/3 < —1 which means we live in a phantom epoch of 
the Universe. 

As a special case, ci = 4 which is a special form discussed [27| we have 



H' = ^-3Hl (12) 

R = 2 ^"^ , I -. (13) 

J% + hR{^-2H^) d r Jj^ + bRi^-2H^)i 

R = {2 ^^ . \ -}^ |2 ^^ ^_\ -] ■ (14) 



w ^R 



Thus ([7]) converts to the following form 

2 J% + bR(§ - 2H^) d J% + bR(§ - 2H^) 

1 2 J% + bR(§ - 2H^) „ R o o 

+-/ = -nsn + fr) + (2 ^^ ^_^^ ^-ffRR + 9(| - 3HI)T^ hr 



We solve equation (|15p using an ansatz f{R,T) = F{T). The solution of (|15|) is 

F(r)=C.T'.+C.T'=.+r-^^^|^^, (16) 



where Ci and C2 are two constants of integration and 



1 -A-9R + 108H'^ + v^l6 + lUR - 1512F^ + 81i?2 . 


- 19URH'^ + 11664^4 


36 R- 9H^ 


-14 + 9R- 108H'^ + VI6 + UiR - 1512H'^ + 81R^ - 


- 19AARH^ + 11664//4 



(17) 



36 R-9HI ^ ' 

If we choose f{R,T) = C3 a pure constant, then Eq. (fT5]) gives C3 = 327rT where T is also 
a constant. Furthermore if f{R,T) = F{R) then Eq. (jlSp cannot be solved analytically or 
numerically. 

D. de Sitter Universe in f{R,T) Gravity 

If we live in a Universe filled by dust and the scale factor increases exponentially with time as 
a{t) = aoe-'^"*, then Hubble parameter is constant. Such model which has been proposed firstly by 



de Sitter, is called the de Sitter Universe. In this solution, it is assumed that the only source of 
matter filling the Universe is dust. Thus it is contained in our treatment on dust formation. For a 
de Sitter Universe in it's static patch we've 

H = Ho, (19) 

R = l2Hl, (20) 

R = R = 0, (21) 

Eq. ([Tj) reduces to 

^RfR + ^f = ^T(87r + /t) + mlT^fTT " 3r(| - S^^) h- (22) 

Solution of equation (p2]) is 



327rr 

'0 
where Fi and F2 are two arbitrary constants and 



f{R, T) = FiT"'' + FaT'^^ - jf^, (23) 



*-^(-^^#^)^ 



.^^-.(WglM). p. 



E. Einstein Static Universe in f{R,T) Model 

In this case since H = 0, thus we have H = R = R = R = 0. Thus we have 

i/ = |T(8^ + /T). (26) 

The solution for this purely T dependence equation is 

f(^T) = -^ + C,T^/\ (27) 



IV. MODELS FOR p = u;p 

Starting from (|4]) and assuming EoS p = up, we have 

Qfiu = -2(1 + uj)pu^Uy + p.g^y (28) 



The trace of (|28j) is given by 

G = 2p{uj - 1). (29) 

Since T = p{l — 3a;) thus, we have 

Now we rewrite (j3]) in the following form 

Rfn + 3n/fi - 2/ = SvrT - T/t - ^^^^/t, (31) 

which can be written in suitable form as 

RfR + SD/fi - 2/ = 8^r - u'TfT, (32) 

where w' = —yz^- We will search for exact solutions of f{R,T) in the following cases. 

A. Solutions in the form f{R, T)^ R + 2/(T) 

Substituting this form f{R,T) = R + 2/(T) in ([32]) we obtain 

-2R- Af{T) = 8ttT - 2u'TfT. (33) 

There is only one possibility to get an exact solution: when R is constant. In this case, the solution 
for (p3]l reads 

f{T) = ciT^/^ + ^ 4(j_2) • (34) 



Thus we have 



/(.,r,....(..v.,?e_^3^±i5^)^ 



B. Solutions in the form f{R,T) = fi{R) + J2{T) 



(35) 



Substituting this form of f{R,T) = fi{R) + /2(T) in d^D we obtain 

Rf[ {R) + 3af[{R)- 2/i {R) = SttT - uj'Tf'^ (T) + 2/2 (T) . (36) 

Note that the left and right hand sides of (j36|) are functions of R and T respectively. Thus solving 
the T dependent part of (p6|) . we obtain 
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FIG. 1: Evolution of Hubble parameter (|^^ against redshift z. We choose C\ = 6, iJ(0) — H^ — 74.2, 
^(z = 0) = 22.26. 



The solution for the R dependent part of the (f36]l is complicated. Indeed we must solve the 
following equation 



R![{R) + 3(9ii/((i?) + 3Hdtf[{R)) - 2h{R) = ci, 



(38) 



or the following equivalent form 



Cl 



Rf[{R) + 2.{Rf'l{R)+R^f';'{R) + mRn{R))-2h{R) = -i, 

b 



(39) 



One simple but interesting solution is obtained by taking Ricci curvature to be constant. From 
(l39l) we have 



HR) 



Cl 



(40) 



Thus one of the interesting models is 



1 2 levrT + (u' - 2)ci 



(41) 



For this model, we can study the evolution of the Hubble parameter and the scale factor. For 
this purpose, we assume fi{R) = R^ and rewrite (j39p in terms of the Hubble parameter H(z) (and 
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FIG. 2: Evolution of log(a) ^ against time t. We choose Ci = 6, o(0) = 1, d(0) = 74.2, a(0) = 2.73 x 10^ 
a(0) = 4 X 10^ 



using dt 



1 dz • 

l+zTTlz)- 



and the scale factor a{t): 



18H 



JH 

dz 



/dH\^ 



£H 



dH^ 



6(1 + .)(_) - 6(1 + ^)H^ + Zh{x2H- - 6Hil + .)-) 



£1 
6 ' 



nd'^a fda\'^d'^a fd^ax"^ da d^ ci , 

a^-T-5( — I ^-T + a[—Tr] +3a— ^-^ = ^a^. 



dt^ 



\dtJ dt'^ \dt'^J 



' dt dt^ 



6 



(42) 
(43) 



The two equations (I42p and (I43p are numerically solved and the result is shown in figures 1 and 2. 
The initial conditions are determined using the definitions of Hubble parameter and jerk parameter. 
Specifically a(0) = 1, Hq = 74.2, d(0) = HqUo = 74.2, go = -0.67, jo = -0.98 [35], a(0) = -H^qo, 
a(0) = —JqHq = 4 X 10^ [35]. Note that our model is compatible with the BAO data [36(] for the 
Hubble parameter. 



C. Reconstruction Using Chaplygin Gas 



The Chaplygin gas (CG) has the equation of state p = ^. There are some extensions of this 
model like the generalized Chaplygin gas but we restrict ourselves only to the CG case. Thus we 
have quantities 



e 



= -2T + 4-, 
P 

P 



(44) 
(45) 



10 



Solving (j45p for p and inserting the solution in ()44p and by defining a new parameter as 12^4 = Tq 
we have 



e ^ em = 2 t^-st^-st^^ 

^ ^ 3 T+y^W^T^ 



Using ([16|) in dH), we have 



Rfn + ^afn - 2/ = SvrT - T/^ - |^o ^J' ^J^Z-^fT. 



(46) 



(47) 



3 T+^T^-Ti 

We are searching for constant curvature solutions in which R = Rq ^ D/j? = 0. Thus solution of 
dUD is 



/(T) = Ro + e 



2/ 



dT 
T+e(T) 



C + Stt 



TdT 



g{T){T + Q{T)) 



(48) 



D. /(i?,T) Models for Scalar Field 



We know that there is an important duality between f{R) models and scalar fields 2j]. The 
Lagrangian for a scalar field, which is minimally coupled to the background reads 



1 
Here w is a free parameter. The corresponding expression for stress-energy tensor is 

Now, the expressions for B and T read 



Eliminating the term (p-a'P'"' from (j51l) and (I52p we have 



G = -6T. 



We can rewrite (jl]) using (f53l) 



RfR + Safn - 2/ = SttT + 5r/T. 



(49) 

(50) 

(51) 
(52) 

(53) 

(54) 



Again, we limit ourselves to constant curvature solutions in which R = Rq (choosing f{R,T) — > 
R + f(T)). Then, the general solution for ([54]) is 

f{R,T) = R-^T + -^. (55) 



11 



V. CONCLUSION 

In summary, we reconstruct cosmological models in the framework of a newly proposed model of 
f{R, T) gravity, where R is the Ricci scalar and T is the trace of the stress-energy tensor. We show 
that the dust fluid reproduces ACDM, Einstein static Universe, de Sitter Universe, phantom-non- 
phantom era and the phantom cosmology. Further, we reconstruct different cosmological models 
including, Chaplygin gas, minimally coupled scalar field, with some specific forms of f(R,T). 
We demonstrated that for one case f{R,T) = R^ + /{T), we found that the behavior of Hubble 
parameter H{z) and scale factor is compatible with the observational data of BAO for small redshift 
z <2. 
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